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Abstract 

Using the extended Poincare-Lighthill-Kuo (PLK) reductive perturbation method, which incor- 
porates the phase-shift variations, it is shown that common features on propagation and head-on 
collisions of ion-acoustic waves exist for a magnetized plasmas of different inertial-less particle 
distributions. For instance it is remarked that, the soliton amplitude is always independent of 
magnetic field strength while strictly depends on its angle regarding the propagation direction. 
Both types of solitons (compressive or rarefactive) are shown to exist which are defined through 
the critical angle 7 = tt/2 or other critical values depending on plasma fractional parameters. 
These critical plasma parameter values also define the sign of head-on collision phase shift. Fur- 
thermore, it is proved that for a given set of plasma parameters there is always a relative angle of 
propagation regarding to that of the magnetic-field for which the soliton width is maximum. Cur- 
rent findings applies to a wide range of magnetized plasmas including those containing background 
dust ingredients or two-temperature inertial-less particles and may be used to study laboratory or 
astrophysical magnetoplasmas. 

PACS numbers: 52.30.Ex, 52.35.-g, 52.35.Fp, 52.35.Mw 



1 



I. INTRODUCTION 



Historically, the first experimental evidence of the ion-acoustic solitary excitations in 
plasma has been established by Ikezi et al. in 1970 fl]. However, due to the wide and 
inevitable applicability in fast growing plasma technology and astrophysical sciences, the 
theoretical prediction of nonlinear density excitations in ionized environments has been 
taken place much earlier 2|, |3j. Since then, several methods such a pseudopotential and 
reductive perturbation techniques [4J has been employed to explore interesting properties 
of such waves, in particular, electron-acoustic, ion-acoustic or dust-acoustic nonlinear struc- 
tures in both classical and quantum plasmas with different species and distributions 



More recently, using magnetohydrodynamics (MHD) model, it has been shown that despite 
the differences in species type or equations of state (EoS), essential dynamics of nonlinear 
excitations in plasmas in the presence of an ambient magnetic field possess common fea- 



tures 
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321 ] . During the past decades, however, a new motivation has started towards the 



study of such nonlinear structures in positron containing plasmas due to their presence in 

J L 

Van Allen radiation belts and near the polar cap of fast rotating neutron stars [33j, active 
galactic nuclei 34], quasars and pulsar magnetosphere 35, 36] and solar atmosphere j^. 



It is also strongly believed that the electron-positron plasma excitations may be the key to 
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39]. 



understanding of the evolution of early Universe 

In recent years many attention have been payed to the generalized Boltzmann-Gibbs 
thermodynamics first formulated by Tsallis 40]. It has been proposed that, this generaliza- 
tion of particle distributions, for instance, may arise due to transfer of the plasma particles 
through the strongly turbulent, non-integer dimensional, and irregular media which may be 



encountered in astrophysical situations 



4=11 ] . This anomalous irregularities in space plasmas 



may also be due to the effects of unknown external forces acting in astrophysical environ- 
ments or due to the wave-particle interaction leading to Tsallis-like distribution of plasma 
species causing a high-energy tail to appear in the distribution function of the particles. 
However, the origin of high-energy superthermal (non-Maxwellian) charged particles, ob- 
served in solar wind, magnetosphere, interstellar medium and auroral zone 42J, |43J, is one 
of the unsolved problems in the field of space and astrophysical plasmas. The application 
of Tsallis-like (Kappa) distribution in modeling the space plasmas had been first suggested 
long before the discovery of Tsallis thermodynamics by Vasyliunas in 1968 44J, and was 
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later adopted by many authors in various physical contexts. There has been extensive the- 
oretical work focused on the effects of superthermal particles on different types of linear 



and nonlinear collective processes in plasmas 45M58| . Among other applications of the su 



perthermal q-nonextensive electron distributions are the interpretation of observations in 
the Earth foreshock and solar wind models with coronal electrons first observed in various 
experimental plasmas, such as laser-matter interactions or plasma turbulence 59]. Unfor- 
tunately, however, there is a lack of investigation of such plasma excitations in the presence 
of magnetic field which is inevitable in astrophysical situations. 

In current investigation we attempt to generalize the previous studies on propagation and 
collision of ion-acoustic solitary waves in three-component electron-positron-ion magneto- 
plasmas employed for various electron-positron distributions to find common features ruling 
the dynamics of such nonlinear structures. We further pay special focus on superthermal 
distribution which is most relevant to space plasmas. The article is organized in the follow- 
ing manner. The basic normalized hydrodynamics equations are introduced in section [Til 
Nonlinear evolution and collision parameters are derived in section IHIl Global and specific 
features are presented in Sees. IIVI and |VJ The final remarks are drawn in section IVIl 



II. HYDRODYNAMIC MODEL EQUATIONS 

We consider a three-component uniformly magnetized plasma consisting of inertial pos- 
itive hot species, say ions and inertial-less negative and positive ingredients, say electrons 
and positrons the masses of which is ignored compared to inertial components. Furthermore, 
we assume inertial particle temperature to be much lower compared to that of inertial-less 
ones, i.e., Tj <C ^+,-- The complete set of three-dimensional normalized MHD equations, 
may be written as 

^ + V • (rijUj) = 0, Ui = iui + jvi + kwi, 

^ + (m ■ V)u< + V<p + fVn, + w(u, x k) = 0, (1) 
VV = (n_ -n+ - th), 

where, 7i», and u) = oj ci /uj p i (u Ci = eB /rrii) are the ion density, velocity and plasma 
ion-frequencies (defined below), respectively. Also, the quantity a = Tj/T_ is the fractional 
ion-temperature and n+ are the number densities of inertialless charged ingredients. The 
dimensionless set of equations, Eqs. (pP) has been obtained using a general scaling defined 
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below 

v _ . . 

V — >■ t- V, t —7- — , n — > n-on, Ui — >■ qui, tp — )■ — <p, (2) 
Aj e 

where, Co> P j = A/47re 2 72 /mj, Aj = Ci/oo P i and q = ^J~ujrrii are the characteristic plasma 
frequency, ion gyroradius and sound-speed the values of which along with the parameter 
v will be defined later based on the charge distribution. Note that, we only consider the 
low-frequency ion-acoustic solitary waves (IASWs) in which u P i <C u>d or equivalently when 
the ion thermal velocity is much less than the value u C i/k. The last equation in Eqs. ([1]) 
may be Taylor-expanded using inertial-less charged particle distributions to give 

Acp = (5 - rii + a 1 (a,n)<f + a 2 (a,/i)v? 2 + (3) 

where, 

n_o T+ 

The the equilibrium charge neutrality condition is given by Poisson's relation as 

a + /3 = 1, P = —, (5) 
n_ 

To find stationary localized solutions we may transform the normalized plasma equations 
(Eqs. (CQ)) to the appropriate stretched coordinate which admits the seperation of variables 
permits elimination of secular terms. The transformation which leads to the desired evolution 



63| 



equations and the corresponding collision phase-shifts is as follows |60l - 
£ = e(kx + ly + mz - c^t) + e 2 P (r], r) + e 3 P 1 (£, rj, r) + . . . , 

v = e {kx + ly + mz- c v t) + s 2 Q (£, r) + e 3 Qi(£, rj,r) + . . . , (6) 

T~ — £ t j — — 

where, the functions Pj and Qj (j = 0, 1,2, ...) describe the phase changes in the traveling 
solitary waves. These quantities will be determined later along with the wave evolution 
equations. It is also assumed that the interacting solitons are initially far apart and travel 
towards each other and described by cosine indices (k, I, m) relative to a magnetic field lines 
with angle 7 defined through 

m = cos 7, 

(7) 

k 2 + I 2 + m 2 = 1. 
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Expanding the dependent plasma variables around their equilibrium values through small- 
ness e parameter, a measure of nonlinearity strength [64], we obtain 
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Stretched state equations are presented in appendix A, isolations of which in lowest-order 
in e leads to the following relations 



+ d v ) nf ] + m/3 + d v ) 
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(i) 



0. 



k/3 + <%) (^(^ + ma (<% + < J - pov. 



(i) 



Z/3 (<% + 0„) y>« + Za + 9„) + /fam} 



(i) 



(9) 



c(3 + 0„) } + m/3 + c^) </? w + ma (9 € + d v ) 



(i) 



n 



(i) 



n- = aiv? (1) , 



from which one obtains the following first-order approximations 
«J 1) = -^[^ (1) (^t) + V 1) (^t)]. 

5 = /3 + at. 

Dispersion relation is, thus, given as 

(3m 2 



(10) 



c 2 — am 2 



ax, 



(11) 



and the corresponding phase-speed c is 



'/? + aid 
ax 



cos 7. 



(12) 



5 



The next higher-order step in isolation of e-power lead to the second-order plasma-state 
approximation of the similar form as above 



+ d v ) nf ] + k (<% + d v ) u\ L> + /(<% + d v ) v?> + m («% 



(i) 



c(3 + d n ) uf ] + k(3 (<% + d v ) <p™ + ka (% + d n ) iif' - favf 1 = 



(2) 



,(2) 



d v ) w. 

,(2) 



(2) 



0. 



c(3 + d v ) vf ] + 1(5 + ^) + /a + d v 



(2) 



(2) 



.-,,,(2) 
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c/3 + 0„) u;f ) + m(3 (<% + 9„) <^ + ma (<% + 0„) 
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which immediately yield the second-order components 
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fll [^ 2 )(e,r) + ^ 2 )(7 7 ,r)], 



^ = m \- d ^ {2) to> r )] + § te^ & r ) + 9 ^ {1) (v,t)] 
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Note that the notations ip^(£, t) and ip^(r], r) describe the first-order amplitude evolution 
and r ) an d ¥ {Vi r ) denote the second-order amplitude evolution of distinct solitary 

excitations in the oblique directions r]± and £j_ (j]± = — respectively. For further simplic- 
ity, we will use the notations tp^ and tp^' instead of p^(^, r) and p^\r], t) in forthcoming 
algebra. 

III. DYNAMICS OF HEAD-ON COLLISION AND PROPAGATION 



The third-order approximation for density component is obtained via isolation of terms 
in higher-order e-power. This is accomplished by solving coupled differential equations in 
this order, making use of dispersion relation (Eq. [TTT) and the first- and second-order plasma 
approximations defined in previous section 
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F(e,r) + G(77,r) 



(15) 
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in which, F(£, r) and G(j], r) denote the homogenous solutions of differential equations with 
the other undefined coefficients in Eq. (fT5l) defined as below 



K 
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A = —L TJ= cos i> 



4( / 8+ai<7) ' 
2a 2 B 2 +a\cj 

R _ /3<D 2 +(ai+/3)(/3+ai<T)sin 2 7 

o -2 / 3 7 '' 

iV= x y sec 7 , (16) 

i?i = 2(3,2 — Jr> 

P _ 4ai(/3+aio-) 
-C/2 — ^ > 

C = 0.2 + 

(ai+/3)(/3+aicr) . 2 1 

To fully determine the dynamical characteristics of the interacting solitary waves we elimi- 
nate the secular terms appearing in Eq. ( fl5l) . which leads to two coupled differential equa- 
tions for each wave, namely 

P o(v,r) = § ; [ tpVdq, (18) 



dr v dr) drj 3 



(19) 



Q (£,t) = ^-J <p™d£, (20) 

On the other hand, the single-soliton solutions for Eqs. ( IT?]) and ( 1T9]) require the perturbed 
potential components and their derivatives to vanish at infinity, i.e. 



c hm o {^ 1) ,^,^} = 0, C = e,* (21) 



with following solutions 



- ^so 



(22) 
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(1) frfi 

^ ^ 2 (^)' (23) 

where, <^o and A represent the soliton amplitude and width, respectively, and Uq is the 
Mach-value. The collision phase-shifts of solitary excitations are obtained using Eqs. f|T8|) 
and ([2QD and the KdV solutions (Eqs. and as 

^o(^r) = f^ A,tanh(^), (24) 



Q (£,r) = §^ A f tanh(^). (25) 

The overall phase-shifts is then obtained by comparing the phases of each wave long before 
and after the collision 



AP 



post— collision past— collision 



lim [e(kx + ly + mz + ct)] — lim [e(fca; + ly + mz + ct)], 

5=0,r;-s>+oo £=0, n^-oo ' (26) 

AQo Q post— collision Q past— collision 

lim [e(kx + ly + mz — ct)] — lim [sffcr + + tuz — ct)], 

where, the quantities AP and AQ denote the overall phase-shifts of solitons "si" and "s2". 
Finally, making use of Eqs. ( )24"1) . ( 1251) and (JHJ), leads to the following expressions for the 
overall phase-shift in head-on collision 



AP = -e 



2a 2 /3 2 -afo- 
4ai/3(/3+aicr) 
2a 2/ 3 2 -afcr 



IV. GLOBAL FEATURES 

It is clear that the three-component model used here reduces to the simple two-component 
(electron-ion) case for — 1 and to the cold-ion case for a = 0. It is also remarked from 
KdV coefficients, Eqs. ffTB"j) . that in general the solitary wave amplitude only depends on the 
angle of the ambient field and is independent of its strength while the soliton width depends 
on both parameters. Close inspection reveals that the soliton width reaches a maximum 
value as the field angle varies in the range < 7 < tt/2. The value of the maximum soliton 
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width in this range is given by the following analytical expression 
7 m = arccos 




(/? + ai)(/3 + a x a) 



and, at the cold-ion limit we obtain 



7n 



arccos 



; 1+ 



P + ai 



(28) 



(29) 



It is easily observed that, in general, the corresponding value of 7 m can not exceed the lim- 
iting value of 7 ~ 54.74°. On the other hand, evaluation of the collision phase-shift given in 
Eq. (127|) reveals that, when the ions are warm (a ^ 0) for some set of critical parameters 
the phase-shift vanishes and changes the sign. For instance, the corresponding critical frac- 
tional ion-temperature for which the phase-shift vanishes is given by a cr = 2a2/3 2 / a i- ^ is 
interesting that these critical values are completely independent of the strength of magnetic 
field, uj or its angle relative to the direction of wave propagation, 7. 



V. SPECIAL CASES 



The preceding arguments are valid for all magnetized three-component plasma with one 
positive inertial and two inertial-less ingredients with opposite charges and can be evaluated 
if the coefficients a\ and ci2 are known. It can be remarked that the general calculation used 
above may be easily extended to magnetized plasmas containing extra immobile background 
species with either charges. However, in order to investigate the effect of various plasma 
parameters on propagation and head-on collision of ion-acoustic solitary waves, we presented 
in subsections below some specific cases which might be of interest. We also emphasize on 
the Kappa (Tsallis) distribution which might be of great importance in astrophysics. 



A. Maxwell-Boltzmann Distribution 



To start with we consider the most simple case of Boltzmann electron/positron. For 
Maxwell-Boltzmann electron/positron distribution in classical electron-positron-ion plasma, 
we have n e = and n p — (1 — ^e^^ with a = Ti/T e , fi = T p /T e and v = k^T e . Then we 
derive a± = 1 + /u(l — (3) and 02 = (1 — yU 2 (l — /3))/2. Characteristics of solitary ion-acoustic 



waves in such plasma has been considered in Ref. 65j which is consistent with the general 
features mentioned in previous section. 



B. Zero- Temperature Thomas-Fermi Distribution 



For the Thomas- Fermi electron/positron distribution in completely degenerate electron- 
positron-ion plasma, we have n e = (1 + <£>) 3//2 and n p — (1 — 0)(1 + /i<£>) 3//2 with a = Ti/T Fe , 
fx = T Fp /T Fe = (1 — /3) _2//3 for nonrelativistic distribution and fi = T Fp /T Fe = (1 — /3) -1 ^ 3 
for ultrarelativistic electron/positron distribution and v = kBT Fe . Then we derive a\ = 
3(1 + /i(l — /3))/2 and <i2 = 3(1 — fi 2 (l — /3))/8. It is observed that, in this case the value 
of a cr (defined above) is 1/9-times that of the Maxwell-Boltzmann electron-positron-ion 
plasma. Dynamics of solitary ion-acoustic waves and their head-on collision in this plasma 
with both ultrarelativistic and nonrelativistic electrons and positrons has been considered 



in Ref. 



661 ] and agrees with the general features mentioned above. Also, similar features 



can be observed for Fermi-Dirac electron-positron-ion plasma considered in Ref. 



67]. 



C. Partially Degenerate Thomas-Fermi Distribution 

For the Thomas- Fermi electron/positron distribution in partially degenerate electron-ion 
plasma, we have n e = (1 + if) 3 / 2 + T 2 (l + p)- 1 / 2 and = 1/(1 + T 2 ) with a = Ti/T e , /i = 
and v = k B T e . Then we derive a± = (3 + T 2 )/2 and a 2 = 3(1 — T 2 )/8. It is observed that, 
in this case the value of a cr solely depends on the electron temperature, T e . 



D. Tsallis Distribution 



Tsallis [40] in 1988 has proposed a generalization of Boltzmann-Gibbs statis- 
tics based on multi-fractal concept of probability. The observation based 
on interstellar plasma velocity distribution measurments confirms that non- 
Maxwellian distributions are common in the solar wind and in the planetary 
magnetospheres where the velocity distribution have a Tsallis-like power-law 
tail at high energies {681 ]. The Tsallis velocity distribution is a convenient ex- 
tension of the well-known Boltzmann velocity distribution in the sense that 
it reduces to the ordinary Boltzmann distribution in the limiting case of the 



spectral index [69J]. However, there are various equivalent versions to character- 
ize the q-nonextensive (Tsallis) velocity distributions usually emplo yed in the 
literature. Here, we consider the distribution adopted in the Refs. [70-72] for 
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inertialless ingredients keeping non-Tsallis distribution for ions. The reason is to 
avoid the Landau damping by keeping the fractional ion to inertialless-particle 
temperature o much smaller than unity. Hence, for the superthermal elec- 
trons/positrons in electron-positron-ion plasma, we have for the spectral index 



values q > 1 [73J| n e = [1 + {q - i)^](<H-i)/2(«/-i) and n p = a[l + {q - l)/^]^ 1 )/ 2 ^ 1 ) with 
a = Ti/T e , fi = T p /T e , a = 1—/3 and v = ksTg. Hence, we derive a\ = (l + g)(l + a/z)/2 
and a 2 = (3 — g)(l + q)(l + a/i 2 )/8. It can be easily verified that the special case 
of q = 1 corresponds to the well-known Maxwell-Boltzmann distribution given 
in Sec. ??. Let us now consider in detail the dynamics of propagation and 
head-on collision of general Tsallis-distributed magnetized plasma which may 
be encountered in various astrophysical situations. 

Figure 1 present the variations of the soliton amplitude with respect to frac- 
tional positron to electron temperature \l for different fractional ion number- 
density f3. It is noted from Fig. 1(a) that for some set of values of //, j3 and 
a the soliton changes from compressive to rarefactive and viceversa. It is ob- 
served from Fig. 1(b) that a critical angle 7 = n/2 and critical fractional positron 
temperature exist which define the shape (brightness/darkness) of the solitons. 

On the other hand, Fig. 2 presents the variations in the soliton width with 
respect to different plasma parameters. It is clearly remarked from Fig 2(a)that 
the soliton width decreases with increase of fi while it increases with increase 
of the fractional ion temperature, a. Also, Fig. 2(b) indicates that the soliton 
width has a maximum value in the range < 7 < 71/ 2 regarding the angle of 
the propagation with respect to that of the field, in agreement with the general 
features mentioned in previous section. It is further remarked from Fig. 2(b) 
that the increase in the value of spectral index, d leads to decrease in the soliton 
width for angles of the ambient magnetic field. 

The head-on collision phase-shift and its variations with respect to various 
plasma parameters are given in Fig. 3. It is clearly remarked from Fig. 3(a) 
that the sign of the collision phase-shift can be positive or negative, in general, 
depending on the chosen plasma parameter set. A positive phase-shift indicates 
that the post-collision parts of the soliton moves ahead of the initial trajectory, 
whereas, a negative phase-shift denotes that the post-collision parts of the soli- 
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ton lags behind the initial trajectory [74]. It is also confirmed from Fig. 3(b)that 
the spectral index has significant effect on the value and the sign of the collision 
phase-shift. Figures 3(b) and 3(d) reveal that the sign and the value of the colli- 
sion phase-shift with varied strength of the ambient field can be much different 
below and above a critical plasma parameter, let say in this case, (3 cr . 



VI. SUMMARY 



Using an extended multiple scales technique, which includes the phase-shift variations, 
we showed that in a magnetized plasmas with diverse inertial-less ingredient distributions 
common rules apply on propagation and head-on collisions of ion-acoustic waves. In general, 
the soliton amplitude is independent of magnetic field strength but strictly depends on its 
angle regarding the propagation direction. It was also shown that the type (dark or bright) 
of the solitons are defined through the critical angle 7 = 7r/2. Moreover, it was shown that 
one or more critical plasma parameter values may exist defining the sign of collision phase 
shift. Current study may be applied to a wide variety of magnetized plasmas including those 
containing background dust ingredients or two-temperature inertial-less particles. 



Appendix A: Stretched Plasma Equations 



(Al) 



P 2 drij _ dni_ _ 2„ dQp dn t , drn _i_ r 2 r 8Pp dm , j drijUj , 
C dr u <9£ c u d£, drj " r u &q c u dri <9£ " r ^ 9£ " r 

r 2i„ dQo drijUj , j dnjUj , 2i„ dPg dnjUj i ] dn0H , ^l dQo dnjVj , 
fc h d£ dri " r h dri " r fc h &q d(, " r 1 9£ " r fc 1 9£ St? _r 

i dTHVj _|_ 2^ dP dnjVj i m §ivwi , ^l^ dQp dn^Wj , dnjWj i 
9»? 9»y 9f 9£ 9£ drj df) 

r 2 duj _ (hH _ r-2„ dQo duj , dv± _i_ r 2 r dPp &Uj . t, duj , 
fc 9r c 9£ fc ° 9£ 9r; " r ° 9r; ~r fc ° 9r; 9£ ~r ft '" 1 9£ " r 

r 2r, 9Q I U„. duj , c 2u dPp§JM _i_ duj i g 2; , 9Qp 9«» i 

t d£ 9^ ' 977 "t" fc a,, a 5 t ^ -r t ^ ar) -r 

^ t + t + m ^ t + If + mw > t+ (A2) 

p 2™„.,.9_Fb 9^. 4- 4- c-2l9Qo 9^ i i P 2 t, 9P 9y i 

£ mw 4 ^ a? 9^ 9£ 9r, + K 9»; + 5 K 9t? 9£ + 

Ua_drH , 2h,a_dQpdrH _i_ h,o_drH , 2; (7 9Pp 9nj 
rii a£ Hi a?7 or) m or\ oc, 



±f + ... = o, 
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r 2 dvj _ d%H _ r 2 r dQ dvj , dvj_ , z 



+ 



fc AiUj ^ ^ f AiUj ^ -t- fc AtUj a? f tWj ^ -t- fc 9 £ ^ i- 



drj 



dr/ 9§ 



9? 



9£ 9r; 



<9?7 



r 2 rnin ^S-^i- -I- -I- r 2/9Qo 9<£ _i_ /£<£ i ^2] dPp dip , 

t 9j? 9? -r <■ 9? -r t <■ a? 9r) -r t a?7 -r t <■ ^ a ^ t 

^ cr drij _|_ ^.2^ <t 9Qo 9wj , ^a_drti , r-2^ it 8Pq drij , 
n, 9§ n i 9r) rii dr\ 9£ 
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0, 



(A3) 
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FIGURE CAPTIONS 



Figure 1 

(Color online) The variation of soliton amplitude with respect to various fractional plasma 
parameters, such as fractional positron to electron temperature \i for different fractional 
plasma parameters such as the relative ion-temperature a, the relative ion number-density 
f3, the spectral-index q and the ambient magnetic-field angle 7, while the other parameters 
are kept fixed. The values of e — 0.1 and u^ = u V:0 = 0.1 are used for all plots in this figure. 
The dash sizes in all plots are appropriately related to the values of varied parameter. 

Figure 2 

(Color online) The variation of soliton width with respect to various fractional plasma 
parameters, such as fractional positron to electron temperature \i for different fractional 
plasma parameters such as the relative ion-temperature a, the relative ion number-density 
j3, the normalized magnetic field strength u and the ambient magnetic-field angle 7, while 
the other parameters are kept fixed. The values of e — 0.1 and = u v fl = 0.1 are used 
for all plots in this figure. The dash sizes in all plots are appropriately related to the values 
of varied parameter. 

Figure 3 

(Color online) The variation of collision phase-shift with respect to various fractional 
plasma parameters, such as fractional positron to electron temperature /x for different frac- 
tional plasma parameters such as the relative ion-temperature a, the relative ion number- 
density (3, the normalized magnetic field strength Co, the spectral-index q and the ambient 
magnetic-field angle 7, while the other parameters are kept fixed. The values of e = 0.1 
and w^o = u Vt0 = 0.1 are used for all plots in this figure. The dash sizes in all plots are 
appropriately related to the values of varied parameter. 
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